function transect_map(ctr,All);
transect=All.transect;
peak=All.peak;

%Load and plot grid data
R=region_data_all(ctr);
[lon, lat, z]=grdread2(['grdfiles/',R.name,'.grd']);
pl=find(lon>=R.cull_box(1) & lon<=R.cull_box(2));
lon=lon(pl); z=z(:,pl);
pl=find(lat>=R.cull_box(3) & lat<=R.cull_box(4));
lat=lat(pl); z=z(pl,:);

%load region data
fname=['../../generate_profiles/data/data_',R.name,'.mat']; 
load(fname,'C','P');

%plot map
h=surf(lon,lat,z,'FaceLighting','flat');
colormap(jet);
axis(R.cull_box);
if isempty(R.clim); R.clim=[min(z(:)) max(z(:))]; end;
caxis(R.clim);
shading flat;
for ct=1:3,
    light('Position',[-1 1 0]);
end;
h.AmbientStrength = 0.7;
h.DiffuseStrength = 0.7;
h=colorbar; h.Label.String='depth (meters)'; h.Label.FontSize=15; h.FontSize=15; 
font(gca,15); 
h=xlabel('longitude'); font(h,16);
h=ylabel('latitude'); font(h,16);
box on;

%plot transects
h4=axis;
ca=2;
for ct1=1:length(C.T) 
    text(C.T(ct1).ridge.rxy(1)+diff(h4(1:2))/50,C.T(ct1).ridge.rxy(2)-diff(h4(3:4))/50,num2str(ct1),'fontsize',20,'color','k');  
    for ct2=1:2,  %plot one side and then the other
            if ct2==1,
                    %identify and mark BM, and find data younger than this
                    ipl=P.T(ct1).BM(1).ipl;
                    [~,p(1)]=min(abs(C.T(ct1).ilon-P.T(ct1).ilon(ipl)));
                    [~,p(2)]=min(abs(C.T(ct1).ilon-C.T(ct1).ridge.rxy(1)));
                    p=sort(p,'ascend');
                    plbm=p(1):p(2);
                    plot(P.T(ct1).BM(1).lon,P.T(ct1).BM(1).lat,'m+','markersize',25,'LineWidth',4);
            else,
                    %identify and mark BM, and find data younger than this
                    if ~isempty(P.T(ct1).BM(2).idist)
                    ipl=P.T(ct1).BM(2).ipl;
                    end 
                    [~,p(1)]=min(abs(C.T(ct1).ilon-P.T(ct1).ilon(ipl)));
                    [~,p(2)]=min(abs(C.T(ct1).ilon-C.T(ct1).ridge.rxy(1)));
                    p=sort(p,'ascend');
                    plbm=p(1):p(2);
                    plot(P.T(ct1).BM(2).lon,P.T(ct1).BM(2).lat,'m+','markersize',25,'LineWidth',4);
            end;   
        plot(C.T(ct1).ilon(plbm),C.T(ct1).ilat(plbm),'k','linewidth',4); %plot line included in this study in black
        plot(C.T(ct1).ridge.rxy(1),C.T(ct1).ridge.rxy(2),'k.','markersize',36); %mark center point in black
    end;
    
end;
